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REPORT 1031 


A STUDY OF THE USE OF EXPERIMENTAL STABILITY DERIVATIVES IN THE CALCULATION 
OF THE LATERAL DISTURBED MOTIONS OF A SWEPT- WING AIRPLANE AND 

COMPARISON WITH FLIGHT RESULTS 1 

Bv John D. Bird and Byron M. Jaquet 


SUMMARY 

An investigation was made to determine the accuracy with 
which the lateral flight motions oj a swept-wing airplane could he 
predicted from experimental stability derivatives determined in 
the 6 -foot-diameter rolling-flow test section and 6- by 6 -foot 
curved-flow test section of the Langley stability tunnel. In 
addition determination of the significance of including the 
nonlinear aerodynamic effects of sideslip in the calculations of the 
motions was desired. All experimental aerodynamic data 
necessary for prediction of the lateral flight motions are pre- 
sented along with a number of comparisons between flight and 
calculated motions caused by rudder and aileron disturbances. 

In general, the agreement between the calculated and measured 
motions of the airplane considered was good when the effects of 
all control movements were taken into account. The greatest 
disagreement occurred at lift coefficients where Reynolds number 
effects on the experimental derivatives would be expected to be 
high, which for the case considered was for lift coefficients above 
about 0.8 when wing slots were used. The nonlinear effects of 
sideslip for this airplane were not very significant for the motions 
considered, which generally involved sideslip angles less than 10°. 

INTRODUCTION 

For the past few years, numerous investigations have been 
made in the Langley stability tunnel to determine the effects 
of geometric variables on the static-, rolling-, yawing-, and 
pitching-stability derivatives of various airplane configura- 
tions. (See references 1 to 4.) In the past, however, none of 
the experimental data have been compared with data obtained 
in flight to determine its relative worth. The. purpose of the 
present report is to determine the applicability of the 
experimental stability derivatives to the prediction of the 
lateral disturbed motions of an airplane in flight. The 
equations used for calculating the motions are given in the 
appendix. 

A /(-scale model of a swept-wing version of a conventional 
fighter airplane, which was selected because of the large 
amount of flight data available (see references 5 and 6), was 
tested in the 6-foot-diameter rolling-flow test section and 6- 
by 6-foot curved-flow test section of the Langley stability 
tunnel to determine all the stability derivatives which are 


usually considered necessary to calculate the lateral motions 
arising from a disturbance caused by the rudder or the 
ailerons. Comparisons have been made between the flight 
and calculated lateral motions for a wide range of conditions 
in gliding flight. 

A few calculations have been made to determine the 
effects of nonlinear variations of the aerodynamic forces and 
moments with the angle of sideslip. 

COEFFICIENTS AND SYMBOLS 

The results of the wind-tunnel tests are presented as 
standard NACA coefficients of forces and moments. 
Moment coefficients are referred to a center of gravity located 
at 2] .8 percent of the mean aerodynamic chord. The wind- 
tunnel data and motion calculations are referred to the 
stability axes, which are a system of axes having their origin 
at the center of gravity and in which the Z-axis is in the plane 
of symmetry and perpendicular to the relative wind, the 
A"-axis is in the plane of symmetry and perpendicular to the 
Z-axis, and the F-axis is perpendicular to the plane of sym- 
metry. The positive directions of the stability axes and of 
angular displacements of the airplane and control surfaces 
are shown in figure 1 . 


Lift 



Figure 1. — System of stability axes. Arrows indicate positive forces, moments, and angular 

displacements. 


1 Supersedes NACA TN 2013, “A Study of tbe Use of Experimental Stability Derivatives in the Calculation of the Lateral Disturbed Motions of a Swept-Wing Airplane and Compar- 
ison with Flight Results” by John D. Bird and Byron M. Jaquet, 1960. 
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The coefficients and symbols are defined as follows: 
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lift coefficient (Lift jqS) 
maximum lift coefficient 
longitudinal-force coefficient (XfqS) 
lateral-force coefficient ( YjqS ) 
rolling-moment coefficient ( L/qSb ) 
yawing-moment coefficient ( N/qS b) 


moment of inertia about longitudinal principal 
axis 

moment of inertia about spanwise principal 
axis 

moment of inertia about normal principal 
axis 

longitudinal force along A'-axis 
lateral force along Y-axis 
normal force along Z-axis (Lif t = — Z) 
rolling-moment about A-axis 
yawing-moment about Z-axis 
wing-tip helix angle, radians 
yawing-velocity parameter, radian measure 
rolling angular velocity about A-axis 
yawing angular velocity about Z-axis 
linear velocity of airplane along Y-axis 
free-stream velocity along A-axis 
calibrated airspeed, based on sea-level density 
of air 

angle of sideslip; /3= — \p in wind-tunnel 
tests ^tan" 1 

angle of attack of wing root chord line 


a 

* 

it 


S 

A 

2 


angle of attack of thrust line (a R — 1.2°) 
angle of yaw, degrees 

angle of incidence of stabilizer with respect 
to thrust line, positive when trailing edge 
is down 

control-surface deflection, measured in a plane 
perpendicular to hinge axis 
angle of sweepback, degrees 


free-stream dynamic pressure 



S' 

b 

A 

P 

T 

T m 

P 
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Subscripts: 

a 


r 

j 

v 


wing area 
wing span 
aspect ratio ( b 2 /S ) 
mass density of ait- 
time 

time to damp to half amplitude 

period 

tail length 

aileron 

rudder 

flap 

vertical tail 


WIND-TUNNEL TESTS 
APPARATUS AND MODEL 

The experimental static-lateral-stability derivatives, 
rolling-stability derivatives, and yawing-stability derivatives 
were determined from tests conducted in the Langley 
stability tunnel in which rolling or curved flight is simulated 
by rolling or curving the air stream about a rigidly mounted 
model. 

The tests were made on a conventional six-component 
balance system with the model mounted at the flight center 
of gravity which is at 21.8 percent of the mean aerodynamic 
chord of the wing. 

The full-scale airplane (a swept-wing version of a conven- 
tional fighter) had the quarter-chord line of the wing, just 
outboard of the intake ducts, swept back 35°. Some of the 
pertinent airplane characteristics are given in table I. 
Afore details of the airplane may be obtained from refer- 
ences 5 and 6. 

The i-scale model shown in figure 2 and in the photo- 
graphs of figures 3 to 6 was constructed of laminated mahog- 
any, finished in clear lacquer, and all surfaces were highly 
polished. The model propeller had three metal blades set 
at an angle of 28° at the 0.75 radius. All propeller-on tests 
were made with windmilling propeller. The model wing had 
a removable leading edge so that slats of 0 percent, 40 per- 
cent, and 80 percent of the swept span could be used inter- 
changeably to simulate those of the full-scale airplane. The 
top surfaces of the slats were cast to the contour of the air- 
foil and the slats were extended by means of metal brackets 
which also act as fences to reduce spanwise flow along the 
slot. A cross section through the slot and slat is shown in 
reference 5. 
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Figure 2,— Geometric characteristics of ^6-scale model of test airplane. All dimensions 

are in inches. 



Figure 3. — Side view of H-seale model mounted in 6-foot-diumctcr rolling-flow test section of 
Langley stability tunnel. 



Figure 4.— Rear view of H-scale model mounted in 6-foot-diameter rolling-flow test section of 
Langley stability tunnel. 



Figure 5. — Close-up of 40-percent leading-edge slots on H-scale model. 



Figure 6.— Close-up of 80-perccnt leading-edge slots on H-scale model. 
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The wing had a plain trailing-edge flap with a chord of 
15.1 percent of the wing chord measured perpendicular to 
the hinge axis. The gap was sealed for all tests. As in the 
case of the airplane the main wheels of the model were fixed 
for all tests; whereas the nose wheel and nose-wheel doors 
were removed for all flaps-up tests. 

Shown in figure 2 are the two ventral fins tested on the 
model. The large ventral fin was used for all tests except a 
few with the 80-percent-span slot configuration for which 
the small ventral fin was used. 

TEST AND TEST CONDITIONS 

Trim tests.- — Model trim lift coefficients of 0.33, 0.55, 0.76, 
and 0.95 were selected as representative of those obtained in 
flight tests. The angle of incidence of the horizontal tail 
was measured with respect to the thrust line. 

In order to determine the trim angles of the horizontal tail 
for the previously mentioned lift coefficients, tests were 
made through the angle-of-attack range with the horizontal 
tail set at —5°, —3°, and 1° incidence. From these tests 
the trim angles of the horizontal tail were determined. (See 
table II.) 

Static tests. — In order to determine the static-stability 
derivatives Ci , C„ , and C Y , the model was tested at 
■/-= ±5° through an angle-of-attack range of «=— 2° to 
a— 23° for the flaps up (trim C L —Q. 33) and a= — 2° to 
a=18° for the flaps down (trim C L — 0.7C) for each of the 
slot configurations. Tests were also made at all selected 
test trim lift coefficients through an angle-of-yaw range of 
ip= ±20° to determine the variations of C ,, C„, and C y 
with for all slot configurations. 

Tare tests were made for the 40-percent-span slot configu- 
ration (flaps up and down). The effect of the slots on the 
tares was assumed to be small; therefore, the tares for the 
40-percent-span slot configuration were applied to all con- 
figurations. The Mach and Reynolds numbers for the tests 
were 0.17 and 1.01 XI 0 6 , respectively. 

Rolling-flow tests. — Tests were conducted in the 6-foot- 
diameter rolling-flow test section of the Langley stability 
tunnel, wherein rolling flight of the model was simulated by- 
rotating the air stream. The model was mounted rigidly 
on a conventional support strut. Details of this test pro- 
cedure are given in reference 1 . 

All slot configurations were tested through the angle-of- 
attack range with the flaps up (trim C L =0.33) and with the 
flaps down (trim C L =0. 76) at helix angles pb/ 2V of 0, 
±0.0253, and ±0.0757 radian. The slopes of C u C„, and 
Cy plotted against pb/2V are the derivatives C, , C„ , and 
CV . The 40-percent-span slot configuration (flaps up and 
down) was tested at the selected trim lift coefficients for the 
previously mentioned values of pb/2V from ^=0° to ^=20° 
to determine the variation of C t , C„ , and C y with }/. The 

V V V 

Mach number and Reynolds number for the rolling-flow 
tests were 0.17 and 1.01 X10 6 , respectively. 

Yawing-flow tests. — Yawing-flow tests were conducted in 
the 6- by 6-foot curved-flow test section of the Langley 
stability tunnel. In this section, curved flight is simulated 


approximately by directing the air in a curved path about 
a fixed model. 

All slot configurations were tested in curved flow through 
the angle-of-attack range with the flaps up (trim Cl— 0.33) 
and with the flaps down (trim ^=0.76) at values at rb/2V 
of 0, —0.039, — 0.082, and —0.108. The slopes of C t , C n , 
and C Y plotted against rh /2 V are the derivatives C t , C n , 
and Cy t - 

The 40-percent-span slot configuration (flaps up and down) 
was tested through an angle-of-yaw range of ^=±20° for 
values of rb/2V of 0, —0.039, —0.082, and —0.108 at the 
previously mentioned trim lift coefficients to determine the 
variation of Ci r , C„ r , and CV r with The values presented 
herein are the average of the results at corresponding posi- 
tive and negative angles of yaw. 

The 40-percent-span slot configuration (flaps up) was 
tested at a trim lift coefficient of 0.33 through the angle-of- 
yaw range with the propeller off. 

The yawing-flow tests were made at a Mach number of 
0.13 and a Reynolds number of 0.8 X10 6 . 

CORRECTIONS 

Approximate jet-boundary corrections based on methods 
derived for unswept wings were applied to the angle of 
attack, longitudinal-force coefficient, and rolling-moment 
coefficient, and a blocking correction of 1.01 was applied to 
the dynamic pressure. 

Corrections for the effect of the support strut have been 
applied to C x , C L , C h C n , Cy, C , + , C^, and C Y+ . In rolling 
flow and curved flow, accurate tares were difficult to obtain; 
and, as a result, the derivatives C\ , C„ , (\- , C\ , C„ , and 

P P P r r 

C y are not corrected for the effects of the support strut. 

The derivative (\ was corrected for the effective pitching 
velocity, which exists when the model is tested at an angle 
of yaw, by the following equation: 

Ci —Ci ’ cos iA±/(±L,A,^) 

V V 

where (\ is measured about the wind axis, and f(A,A,\p) 
which is small as compared with C L ' cos \p, is a function 
determined by use of the methods of reference 4. Corre- 
sponding effective pitching corrections were not applied to 
the derivatives C„ and C Y . 

p V 

A correction was also applied to the derivative C Y to 
account for the error caused by the cross-tunnel static- 
pressure gradient which is associated with curved flow. 

EXPERIMENTAL RESULTS 

The experimental data are discussed briefly with reference 
to the effects of the slots and angle of yaw on the aerody- 
namic characteristics of the model, because the effects of 
these variables on the rotary derivatives have not been 
investigated extensively to date. The figures which present 
the results obtained in the present investigation are listed in 
table III. 
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The basic lift and longitudinal-force data of figures 7 and 
8 are generally in good agreement with larger scale tests of 
another model of the same airplane, as given in reference 7. 

I he main effect of the slots is to extend the linear range of 
those stability derivatives which are largely contributed by 
the wing to higher lift coefficients in a manner similar to the 
effect of slots on the lift curve, of a wing. ..... (See, figs. 7 to 10 

and 12, 13, 15, and 16.) One significant effect of the slots 
is on the damping in yaw <7, v which increased as the slot 
span is increased. When the 80-percent-span slots are used 
Cn r is increased (over that of the unslotted configuration) 
by about 25 percent at C, = {). (See fig. 15 (a).) When the 
flaps are deflected (fig. 15 (b)), the effects of the slots on the 
yawing-stability derivatives are not as great as when the 
flaps are retracted. 

A comparison of figures 11(a) to 11(f) shows that C n 
varies to some extent through the yaw range, CV is approx- 
imately constant, and for any given lift coefficient C , is 


approximately constant between ^= — 10° and ^=10°. As 
the angle of yaw is increased (fig. 14 (a)), C Y tends to de- 
crease, C n remains approximately constant, and C\ de- 
creases slightly. The variations are similar when the flaps 
are deflected. (See fig. 14 (b).) 

The tests for the determination of the variation of G v 

... r r> 

C„ r , and C, with \p were made for negative values of rb/2V 
only. Results for positive values of rb/2V and positive ^ 
were obtained by assuming that the model was essentially 
symmetrical about the XZ-plane and by utilizing the results 
for the corresponding opposite angles of yaw and rb/2V with 
regard for signs. This procedure amounted to averaging the 
derivatives for corresponding positive and negative angles 
of yaw. 

In general, as the angle of yaw is increased from 0° to 20°, 
the damping in yaw C n ^ is increased by about 15 percent, and 
CV and Ci r are decreased slightly (fig. 17 (a)). Deflection 
of the flaps or removal of the propeller does not appreciably 
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chang e the variations of the yawing-stability derivatives with 
angle of yaw. (See fig. 17.) 

Because C n is largely dependent on the size of the vertical 
tail it is approximately true throughout the yaw range that 


Substitution of the proper values of the span, C n ^, and the 

tail length l, indicates that the trend of the variation of 
C n with \p shown in figure 17 is reasonable. 



Figure 8— Variation of longitudinal-force coefficient and angle of attack with lift coefficient 
for two slot configurations. Trim Cl =0.33; ^=0°; 5/=0°; R— 1.01X10*. 


The relative constancy of C Yv , C„ v , C Yf , Gi t , and Cn r 
with angle of yaw as indicated by figures 14 and 17 and the 
linearity of the curves of C Y , Ci, and C„ plotted against \p 
for angles of yaw up to approximately 10° (fig. 11) were 
factors which indicated that nonlinearities were not of first- 
order importance for this airplane in the calculation of 
motions involving reasonably small variations in Conse- 
quently, most of the motion calculations neglect the effect of 
\f/ on the stability derivatives. 

The results of unpublished tests of swept wings at Reynolds 
numbers to 8.0X10 6 in the Langley 19-foot pressure tunnel 
indicate that the linear part of the curve of plotted against 
C L is increased by an increase in Reynolds number. The 
curves of C t against C L given herein agree well with those 
obtained in flight (references 5 and 6) and in tests of a 
-^r-scale model (reference 7) except at lift coefficients above 

C L —0.8 where the magnitude of the present test values 
decreases for the slotted-wing configurations. The linear 
parts of the curves of the rolling- and yawing-stability 
derivatives are believed to be extended similarly to higher 
lift coefficients if the Reynolds number is increased. 

MOTION CALCULATION METHODS 

The lateral disturbed motions of the test airplane were 
calculated from the aerodynamic data obtained from the 
tests described in the section entitled “Wind-Tunnel Tests.” 
The mass and dimensional characteristics of the airplane 
are given in table I, and a tabulation of the flight conditions 
for which lateral disturbed motions were calculated for com- 
parison with the flight motions is given in table IV. Most 
of the calculations involved dynamic derivatives which were 
constant for a given lift coefficient as is usually employed in 
the theory of small disturbances used in lateral-stability 
calculations. 

Solutions of the lateral equations of motion, given in the 
appendix, were obtained for unit step disturbances in roll 
and yaw by the method described in reference 8. The 
aileron and rudder deflections during the flight motion under 
investigation were then approximated by a series of step 

functions usually at ^-second intervals. The motion arising 

from the control movements was then calculated by 

m r~ "l 

fin , Prcj 01 i n b* \ I ( ^ m — V ')rudd<r ( ~i~ ( Li — n)aileroii(> B/I )roit I 

W = 0L J 

where B n is the value of the unit solution caused by A( T „=1 
or A77 j= 1 at a time T=kn, and 8 m - n is the fraction of the 
unit disturbance applied by the rudder or aileron at a time 
T=k(m—n). The rudder and aileron effectiveness were 
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obtained from reference 7. This procedure is essentially an 
approximate evaluation of Duhamel’s integral and was con- 
sidered sufficiently accurate for those calculations. Refer- 
ence 8 gives a more exact graphical evaluation of this integral. 
The yawing moment caused by aileron deflection and the 
rolling moment caused by rudder deflection were not con- 
sidered of enough significance to warrant their inclusion in 
these calculations. In some cases, however, these factors 
may be of greater significance. 

In a few cases, unit solutions to the equations of motion 
were obtained by a Laplace transform procedure which has 
been adapted for use with automatic digital computers. 
The results were, of course, identical with those presented 
in this report. 


Calculations of the lateral motions for a few cases employ- 
ing a nonlinear variation of rolling-moment and yawing- 
moment coefficients with angle of sideslip and a variation 
of C„ r with angle of sideslip were carried out by use of the 
Kutta three-eighths rule for solving the lateral equations of 
motion. (See reference 9.) All lateral-motion calculations 
were made on an automatic digital computing machine. 

CALCULATED LATERAL MOTIONS 

GENERAL 

The flight records corresponding to the motions calculated 
for this report showed that the motions resulting from right 
and left control movements were not exactly of opposite 
magnitude. This result indicates that there was some asym- 
metry in the characteristics of the airplane, although the 



Figure 9 —Variation oflatcral static-stability derivatives with lift coefficient for three slot configurations. Large ventral On on; propeller on; ^=0°; /i-1. 01X10®. 


901351—52 2 
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wind-tunnel test results indicated no marked asymmetry in 
the model characteristics. Part of the differences in the 
flight motions to the right and left are undoubtedly due to 
variation hi the thrust conditions. Although the flight tests 
were made under approximately zero thrust, no convenient 
means of obtaining this condition was available, and thus 
this requirement was left to the judgement of the pilot. 
Some small part of the difference in the motions to the left 
and right should be attributable to instrument error as the 
film record frequently contained considerable hash which, 
of course, tended to obscure the actual motion record. The 
differences between the motions to the right and left were 
resolved in the present report by presenting both records 
with the sign of one reversed so that the motions were super- 
imposed. The shaded areas in the lateral-motion plots 
(figs. 18 to 23 and 25 to 35) represent the difference between 
the motions to the right and left. The flight motions corre- 
sponding to the calculated motions were obtained from 
references 5 and 6 and related tests. The figures which 
present the results of the lateral-motion calculations arc 
listed in table III. 



Figure 10.— Variation of lateral static-stability derivatives with lift coefficient for two slot 
configurations. Trim Cl = 0.33; ^=0°; f/=0°; #=1.01X10 8 . 


LATERAL OSCILLATIONS 

Lateral oscillations were initiated by abruptly deflecting 
the rudder of the airplane and returning it to neutral equally 
as rapid a moment later. Some aileron waggle caused by 
the floating tendency of the ailerons occurred and was 
accounted for in the calculations by the procedure given in 
the calculation methods. 

The calculated and flight motions are generally in good 
agreement for all motions calculated (figs. 18 to 23) except 
for the condition at 67—0.977. (See fig. 23.) A comparison 
of the flight and calculated periods and time to damp to 
half amplitude (fig. 24) also indicate fairly good agreement 
except for the motions at lift coefficients of 0.977 and 1.169. 

The values of C, used for the calculation of the lateral 

P 

motions were obtained from the curves of C t plotted against 

(fig. 11) wherever possible rather than the curves of C, 

v 

obtained from tests made at ^=±5° (figs. 9 and 10) because 
this procedure is believed to bo more accurate. Although 
the difference between the two methods is generally small, 
such is not always true; and in some cases the calculated 
rate of damping was found to be appreciably affected by 
the difference in C, 

P 

At the higher lift coefficients, the experimental stability 
derivatives deviate appreciably from their initial trends. 
This tendency previously has been referred to as a Reynolds 
number effect and is probably the cause for the lack of agree- 
ment at high lift coefficients between the flight and calculated 
results. References 5 and 6 which present flight tests for 
this airplane show no reduction in <7^ up to the maximum 
test lift coefficient. The rotary derivatives of the airplane 
presumably behave similarly. Evidence of the deviation of 
the experimental stability derivatives from their true varia- 
tion with lift coefficient is observed in figure 20 which presents 
results for a lateral oscillation occurring at 67=0.759. The 
calculated result indicates an excessive response in sideslip, 
yaw, and roll. 

A few additional calculations of the period and time to 
damp, made for other high lift coefficients in flight, indicate 
increasingly poor agreement between the flight and calcu- 
lated times to damp to one-half amplitude with increasing 
lift coefficient. The relatively good agreement between the 
periods of the flight and calculated motions for all lift co- 
efficients, however, indicates that the experimental values of 
G n are fairly close to the correct values. The period of the 
motion is primarily a function of the directional-stability 
parameter C„. (Sec reference 10.) An extrapolation of 
the curves of the derivatives plotted against lift coefficient — 
which amounted to selecting the value of the derivative just 
previous to the break in the curves occurring at maximum 
lift coefficient — was employed for one case at C L = 1.17 but 
failed to yield a satisfactory result. A linear extrapolation 
of the curves in the region preceding the departure from the 
theoretical or linear trend is expected to bo more satisfactory. 
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Figure 11. — Variation of Cy, C n , and Ci with angle of yaw. Large ventral fin on; propeller on; /?= 1.01X10*. 
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(e) 80-percent-span slots; 5/=0° (f) 80-percent-span slots; 5/=40°. 


Figure 11— Continued. 
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RUDDER KICKS 

Rudder lucks were initiated by abruptly deflecting the 
rudder of the airplane and holding this deflection as the 
airplane responded. The records of the flight motions were 
short, usually covering 5 seconds. As in the case of the 
lateral oscillations discussed previously, some slight aileron 
waggle that occurred was accounted for in the calculated 
motions. The agreement between the flight and calculated 
motions is, in general, quite good (figs. 25 to 31). Agreement 
is best at the lower lift coefficients; however at a lift coeffi- 
cient of 0.919 (fig. 28), good agreement is shown between 
flight and calculated angles of sideslip and rolling velocities. 

The flight tests reported in reference 5 showed that at low 
lift coefficients the airplane rolled in response to a rudder 
kick as if the airplane had a decided negative dihedral 


effect. Figure 31 indicates that this peculiar response was 
largely caused by the slight aileron waggle which occurred 
during the motion. In general, for all configurations investi- 
gated, it was found necessary to account for any slight aileron 
movements in order to predict satisfactorily the lateral 
motions of the airplane. 

The motions at (7^=0.794 (fig. 30) give evidence of de- 
parture of the experimental derivatives from the true varia- 
tion with lift coefficient. This result again is the Reynolds 
number effect previously discussed. The flight and calcu- 
lated periods of the lateral oscillation caused by the rudder 
kicks are in very good agreement which indicates again that 
the values of C n used in the calculations were reasonably 
accurate. 



(g) 40-percent-span slots; 5/=0°; propeller off. 


Figure 11.— Concluded. 


(h) 80-percent-span slots; 5/=0°; small ventral fin on. 




0 .2 A .6 .8 A O /.i 

Lift coefficient ; 

Figure 13. — Variation of rolling-stability derivatives with lift coefficient for two slot 
configurations. ^=0°; 5/=0°; 12=1.01X10*. 
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Figure 15.— Variation of yawing-stability derivatives with lift coefficient for three slot configurations. Large ventral fin on; propeller on; ^=0°; R= 0.8X10 6 . 


AILERON ROLLS 

Aileron rolls were initiated by abruptly deflecting the 
aileron control of the airplane and holding the deflected 
position as the airplane responded. The rudder control was 
held as near neutral as possible. The records of motion 
were necessarily short because of the large angles of bank 
assumed by the airplane after a short period of time. 

Comparison of the flight and calculated aileron rolls indi- 
cates fairly good agreement (figs. 32 to 35) except for the 
high-lift-coefRcient condition without nose slots. The sta- 
bility derivatives for the unslotted configuration show a 
departure at fairly low lift coefficient from the initial trend 
of the variation of the derivatives with lift coefficient. As 
previously mentioned, this result is not obtained at the 
flight Reynolds numbers. The calculations for this high lift 
coefficient were carried out with two sets of stability deriva- 


tives. One set was obtained by a linear extrapolation of the 
curves of the derivatives plotted against lift coefficient and 
the second set, by selecting for the value of the derivative that 
value just previous to the final break in the curves of the 
derivatives plotted against lift coefficient. The second pro- 
cedure was necessary because the flight lift coefficient was 
greater than the maximum experimental lift coefficient. The 
linear extrapolation produces the better result, but it can be 
seen that neither set of derivatives was very satisfactory 
although it would appear logical to believe, in view of the 
Reynolds number effects indicated previously, that the linear 
extrapolation of the derivatives would have given fairly good 
results. The aileron roll for the 80-percent-span slot con- 
figuration with flaps down at (\— 1.169 (fig. 35) shows fair 
agreement between the flight and calculated motions so that 
a beneficial effect of the slots in maintaining an unscparated 
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FifiURE 16. — Variation of yawing-stability derivatives with lift coefficient for two slot 
configurations. Trim Cl=0.33; ^=0°; </=0°; /?=0.8X10*. 


flow over the wing to a fairly high lift coefficient is indicated. 

The effects of elasticity of the wing were not considered in 
the calculations of the aileron rolls because of lack of knowl- 
edge of the flexibility of the wing. Rough estimates of the 
effect of elasticity on the rolling velocity for the results of 
figure 33 indicate that the discrepancies shown between flight 
and calculated results would be reduced by about 50 percent 
or more. 

A comparison has been made of the maximum rolling 
velocities obtained by calculation for conditions duplicating 
those in flight, by flight tests, and by calculations for a 
coordinated maneuver in which the yawing velocity and angle 
of sideslip are maintained at zero. For this last case, the 
rolling velocity becomes 


-A Ci 2V 



The results are presented in the following table: 


Figure 

Maximum rolling velocities from 
figures 

| AC, \2V 

Calcu- 

lated 

Flight 

(left) 

Flight 

(right) 

Ci 6 

1 v 1 

32 

0. 62 

0.44 


2. 67 

33 

.77 

.60 

6.53 

1.11 

34 

.78 


.67 

1. 13 

35 

.34 

.33 

.30 

.68 


The comparison indicates that a marked difference may be 
obtained in the maximum rolling velocity by eliminating 
the degrees of freedom of sideslip and yaw. This last 
method of calculation is used frequently in comparing the 
relative merits of various forms of ailerons. 

NONLINEAR CALCULATIONS 

Calculations, in which the curves of (7, and C n against 
were represented not by a single slope but by a series of 
tangents to the curves and in which the variation of Cn r 
with i/' was included, were made for two motions to illustrate 
the effect of these nonlinearities. The results are given in 
figures 36 and 37 along with the flight motion and the cal- 
culations made for linear slopes and constant damping 
derivatives. The effect of the aileron waggle was small 
and was taken to be the same as for the linear calculation. 
These figures indicate that for the angles of sideslip encoun- 
tered in these motions little is to be gained by going to the 
additional effort required to make the calculations for the 
nonlinear case. The nonlinear ’•calculations required ap- 
proximately ten times as long to complete as the linear cal- 
culations. Motions at large angles of sideslip or motions 
for a configuration having more pronounced nonlinearities 
than for the airplane considered herein would undoubtedly 
show the nonlinearities to bo of greater significance. 

SUMMARY OF RESULTS 

The investigation indicated that the lateral disturbed 
motions of the airplane considered herein may be generally 
calculated with good accuracy by the use of experimental- 
scale-model data up to lift coefficients where any appreciable 
effects of the difference in Reynolds number between the 
scale and flight conditions begin to be apparent. This 
effect is usually evidenced by a departure of the scale results 
from a gradual variation with lift coefficient and, for the 
case investigated, corresponds to a lift coefficient of 0.8 for 
the slotted-wing configurations. At higher lift coefficients 
the accuracy in calculating flight motions is progressively 



16 


REPORT 1031 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Figure 17— Variation o( yawing-stability derivatives with angle of yaw. Large ventral fin on; 40-percent-span slots; R=0.8X10«. 


poorer, especially in the time to damp to half amplitude where 
errors as much as 100 percent or more may be encountered. 
For the airplane of this investigation, at these high lift co- 
efficients, the calculated time to damp to half amplitude was 
generally higher than the flight value. In most cases the 
period of the motion could be calculated to within 5 percent 
of that obtained in the flight motion. In order to obtain 
good accuracy in predicting flight motions, proper considera- 
tion of all control movements must be made. 

For most cases, the nonlinear variation of the directional- 
stability parameter, the effective dihedral parameter, and 
the variation of the damping in yaw with angle of yaw could 


be neglected. If, however, the flight motions reach a suffi- 
ciently large amplitude for a long period of time, as in an 
aileron roll for example, greater accuracy in the calculated 
motions should result if the nonlinearities in the variation of 
yawing -moment coefficient, rolling-moment coefficient, and 
damping in yaw with angle of yaw are included in the 
calculations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 1, 1949. 
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per hour; engine idling. 



Figure 19. Comparison of flight and calculated lateral motions resulting from an abrupt 
deflection and release of the rudder. 40-percent-span slots; flaps up; C*,=0.551; Vc— 156 
miles per hour; engine idling. 



Figure 20.— Comparison of flight and calculated lateral motions resulting from an abrupt 
deflection and release of the rudder. 40-percent-span slots; flaps up; Cl =0.759; Vc=13G miles 
per hour; engine idling. 



Figure 21 —Comparison of flight and calculated lateral motions resulting from an abrupt 
deflection and release of the rudder. 40-percent-span slots; flaps down; Cl— 0.524; Vc=160 
miles per hour; engine idling. 
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Figure 22 .— Comparison of flight and calculated lateral motions resulting from an abrupt 
deflection and release of the rudder. 40-percent-span slots; flaps down; Cr=0.801; Vc=128 
miles per hour; engine idling. 



Figure 23.— Comparison of flight and calculated lateral motions resulting from an abrupt 
deflection and release of the rudder. 80-percent-span slots; flaps up; small ventral fin on ; 
Cr=0.977; Vc=120 miles per hour; engine idling. 



Lift coefficient , C L Lift coefficient, C L 

(a) 6 ,= 0°. (b) 5/=40°. 

Figure 24.— Comparison of flight and calculated values of the period and time to damp to half amplitude. 
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Figure 25.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 40-percent-span slots; flaps up; Cfc=0.341; Vc=198 miles per hour; engine idling. 


Figure 26.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 40-percent-span slots; flaps up; Cl=0.556; Vc=100 miles per hour; engine idling. 



Figure 27.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 40-percent-span slots; flaps up; Cl=0.764; Vc* 133 miles per hour; engine idling. 
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Figure 28.— Comparison of flight and calculated lateral motions resulting from a right rudder 
kick. 40-percent-span slots; flaps up; Cl =0.919; T c=120 miles per hour; engine idling. 



Figure 29.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 40-percent-span slots; flaps down; Cl = 0.537; Tc=157 miles per hour; engine idling. 



Figure 30.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 40-percent-span slots; flaps down; Cl = 0.794; Vc=130 miles per hour; engine idling. 



Figure 31.— Comparison of flight and calculated lateral motions resulting from a left rudder 
kick. 80-pcrccnt-span slots; flaps up; Cl=0.278; Vc=228 miles per hour; engine idling. 
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Figure 32. — Comparison of flight and calculated lateral motions resulting from an abrupt 
left aileron deflection. O-percent-span slots; flaps up; C7„=0.983; Fc= 119 miles per hour; 
engine idling. 
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Figure 33. — Comparison of flight and calculated lateral motions resulting from an abrupt left 
aileron deflection. 40-pcrcent-spun slots; flaps up; (7 l=0.598; Kc= 150 miles per houi; 
engine idling. 
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Figure 36.— Effect of nonlinear variation of Cj, C n , and C„ r with angle of yaw on the calcu- 
lated lateral motions resulting from a left rudder kick. 40-percent-span slots; flaps up; 
Cl = 0.764; Vc— 133 miles per hour; engine idling. 


Figure 37— Effect of nonlinear variation of C t , C n , and C„ r with angle of yaw on the calcu- 
lated lateral motions resulting from an abrupt left aileron deflection. 40-percent-span slots; 
flaps up, C!l=0.598; Vc— 150 miles per hour; engine idling. 



LATERAL EQUATIONS OF MOTION 


APPENDIX 


The lateral equations of motions used for most of the cal- 
culations of this report consider the fluctuation of the aero- 
dynamic forces to be directly proportional to the component 
angular and translational velocities associated with the 
axis system as is customary in classical stability theory. 
The modification to these equations to account for the 
nonlinear effects of sideslip is discussed in the text. The 
equations used herein contain the necessary product of 
inertia and control terms and are given as follows: 

Roll 

Ci 6 'S a + Ci lj (i-\--^Ci }i D<p-\-—Ci r D\(/=2ti cos2 V 

+(~°Jsin 2 „] DV + [(^-(^ 0 )] cos v sin „ j 
Yaw 

^„ s S r - , rC„^5 a J rC„ e l3-\--f, lii D<p- s r-f„ r D\ff=2n | ^ y > 'j cos 2 n 

+(^°)'sin 2 ,] D V+[(x°) -(t- 0 ) 2 ] cos * sin * m) 
Sideslip 

Cy s ^S T -\- Cy s ^Sa-\- C y ^-\~2 Cy v D<p-\-— Cy r D\l/-\- C L ip-{- 
C L \j/ tan y = 2n(D(l-\-D\p) 


rate of change of rolling-moment coefficient 
with rudder deflection 

rate of change of rolling-moment coefficient 
with aileron deflection 

rate of change of yawing-moment coefficient 
with rudder deflection 

rate of change of yawing-moment coefficient 
with aileron deflection 

rate of change of lateral-force coefficient with 
rudder deflection 

rate of change of lateral-force coefficient with 
aileron deflection 



<p 

angle of bank, radians 

w 


pSbg 


W 

weight of airplane 

9 

acceleration due to gravity 

k x 

A 0 

radius of gyration about principal longitu- 


° r « 5/3 


dinal axis —yy) 

radius of gyration about principal normal 
axis (V%/-) 

flight path angle, positive for climb 
inclination of principal longitudinal axis of 
inertia with respect to flight path, positive 
when above flight path at nose 
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TABLE I.— PERTINENT AIRPLANE DIMENSIONS AND 
CHARACTERISTICS 


Mass characteristics: 

Normal gross weight, lb. 

Moments of inertia, ft-lb-sec 2 : 



Principal axis (relative to thrust line at nose), deg. 
Center of gravity: 

Location on M.A.C., percent M.A.C 

Relative to thrust line, percent M.A.C. below 


8, 700 


7, 654 
14, 088 
20, 159 
0. 45 up 

21. 8 
14. 8 


Wing: 

Span, ft 33. 6 

Area, sqft 250 

Airfoil section (normal to L. E.) : 

Root Modified 66, 2x-116 (o— 0. 6) 

Tip Modified 66, 2x-216 (a = 0. 6) 

M.A.C., ft 7.79 

Leading edge of M.A.C. (ft behind L. E. root chord) 3. 27 

Aspect ratio 4. 51 

Taper ratio 1. 84:1. 00 

Dihedral, deg 0 

Sweepback (at quarter-chord line), deg ... 35 

Total area, plain sealed wing flaps, scj ft 12. 52 


Ailerons: 

Span (along hinge line, each), ft 8. 75 

Area (rearward of hinge center line, each), sq ft . . 6. 51 


Horizontal tail: 

Total area, sq ft . . . 46. 53 

Vertical tail: 

Fin area (above horizontal fail), sq ft . 13. 47 

Total rudder area, sq ft . ... . ... . 10.26 

Ventral fin area, sq ft: 

Large ventral . 17. 10 

Small ventral (approx.) 8.50 


TABLE II.— TRIM ANGLES JDF^HORIZONTAL TAIL 


Slots (per- 
cent span; 

V 

(deg) 

Ventral 

fin 

Pro- 

peller 

Trim Cj, 

U 

(deg) 

0 

0 

Largo 

On 

0.33 

-0. 10 




.55 

-1. 10 





.70 

-2. 75 

0 

40 

Large 

On 

.33 

.90 




.55 

— . 60 





.70 

-1. 50 





.95 

-2. 25 

40 

f) 

Large 

On 

.33 

-.00 




and 

.55 

-2.70 




off 

.70 

-4. 00 

40 

40 

Large 

On 

.33 

.75 





.55 

— 1.50 





.70 

-3. 20 





.95 

-4.40 

80 

0 

Large 

On 

.33 

0 



and 


.55 

-2. 40 



small 


.76 

-4. 25 

80 

40 

Large 

On 

.33 

.60 




.55 

-1.40 





.70 

-3.00 





.95 

-4.20 


TABLE III— INDEX OF FIGURES 
Wind-tunnel results: 

Static stability characteristics: Figures 

Longitudinal characteristics 7 and 8 

Lateral derivatives 9 and 10 

Variations of C t , C„, and Cy with p 11 

Rolling derivatives: 

Derivatives at p = 0° 12 and 13 

Variations of derivatives with P - 14 

Yawing derivatives: 

Derivatives at ^=0° 15 and 16 

Variation of derivatives with p 17 

Calculated and flight motions: 

Abrupt deflection and release of rudder 18 to 23 

Summary of period and damping values 24 

Rudder kicks 25 to 31 

Aileron rolls 32 to 35 

Nonlinear aerodynamic effects of sideslip 36 and 37 


TABLE IV.— FLIGHT CONDITIONS FOR WHICH MOTIONS 
WERE CALCULATED 


Figure 

Slots 

(percent 

span) 

Flaps 

(dog) 

Ventral 

fin 

i'. 

(mph) 

Cl 


18 

40 

0 

r 

198 

0. 334 


19 

40 

0 


150 

.551 


20 

40 

U 

1 Large 1 

130 

.759 


21 

40 

40 


100 

.524 


22 

40 

40 

1 1 

128 

.80! 


23 

80 

0 

Small 

120 

.977 


25 

40 

0 


198 

.341 


20 

40 

0 


160 

. 556 


27 

40 

0 


133 

.764 


28 

40 

0 


120 

.919 


29 

40 

40 


157 

.537 


30 

•10 

40 

\ Larne i 

130 

. 794 


31 

80 

0 


228 

.278 


32 

0 

0 


119 

.983 


33 

40 

0 


150 

.598 


34 

40 

40 


146 

.598 


35 

80 

40 


110 

1. 169 



Typo of test 


Oscillation 


}Rudder kick 


Aileron roll 
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